To probe quantum criticality with scanning tunneling spectroscopy 
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We investigate the role of quantum coherence in tunneling conductance, where quantum criticality 
turns out to suppress Fano resonance. Based on the nonequilibrium noncrossing approximation, 
we show that the linear tunneling conductance exhibits weak Fano line-shape with sharp cusp at 
zero energy in the multichannel Kondo effect, resulting from incoherence associated with quantum 
criticality of impurity dynamics. In particular, shift of the peak position in the Fano resonance is 
predicted not to occur for the multichannel Kondo effect, distinguished from the Fermi liquid theory 
in the single channel Kondo effect. 
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I. INTRODUCTION 

Recently, scanning tunneling microscopy (STM) has 
been utilized extensively to probe the electronic struc- 
ture of materials with atomic scale spatial resolution. It 
was found that formation of the Kondo resonance gives 
rise to an asymmetric line-shape in the tunneling conduc- 
tance through the STM tip close to a magnetic adatom 
on a metallic surface^— the origin of which is an effect of 
interference between the direct tip-host tunneling and in- 
direct tip-adatom-host one that resembles so called Fano 
resonance X This Fano-Kondo effect has been discussed 
in detail^— A similar Fano-Kondo effect was also inves- 
tigated in the electron transport through a quantum dot 
embedded in a closed Aharonov-Bohm interferometer A 9 - 

Mechanism of Fano resonance implies that quantum 
coherence of impurity dynamics plays an important role 
for the line-shape, where the coherence time scale is esti- 
mated as ~ 1/Tk with the Kondo temperature Tk- An 
interesting question is what happens in the Fano line- 
shape if impurity dynamics becomes incoherent. Such 
a situation is realized in the multichannel Kondo sys- 
tem, where screening of a local moment by conduction 
electrons is overcompensated to drive the local Fermi 
liquid state into a non-Fermi liquid critical state, first 
suggested by Nozieres and Blandin in the multichannel 
Kondo model.— 

Multichannel Kondo impurity systems have been stud- 
ied both experimentally and theoretically with consid- 
erable interests. Recently, the multichannel Kondo 
model was realized artificially in quantum dots^i The 
multichannel Kondo effect was also claimed to occur 
in the quadrupolar Kondo effect^ and in metal point 
contacts. 13 In the theoretical respect the multichannel 
Kondo model has been studied in a variety of controlled 
techniques^ The conformal field theory approach pro- 
vides exact results of the non- Fermi liquid fixed point ^ 
and the noncrossing approximation (NCA), exact in the 
limit of large number of spin flavors and charge channels, 
also gives practically sensible results , 16 ' 17 where univer- 
sal power-law scaling is found in physical responses. Such 
power-law physics distinguishes the critical non-Fermi 



liquid state of the multichannel Kondo impurity from the 
local Fermi liquid state of the single- channel Kondo im- 
purity. 

In this paper we study the Fano-Kondo effect by the 
tunneling current which flows from a single- channel STM 
tip to a multichannel Kondo impurity host. Instead 
of the standard Fano-Kondo resonance in the tunneling 
conductance, one may expect different pronounced fea- 
tures due to interference between the Fano resonance of 
the tunneling current and the overcompensated screen- 
ing of an impurity. Employing the Keldysh nonequilib- 
rium formalism ; 18 i 19 we derive the tunneling current and 
its conductance, where the tunneling current solely de- 
pends on the Green function of an impurity. We cal- 
culate the linear conductance profile analytically at zero 
temperature, based on the nonequilibrium NCA 2 -i2I to 
obtain nonequilibrium Green functions of the impurity. 
A power-law line-shape in the tunneling conductance 
clearly shows the overscreening effect of an impurity. 
Such an effect of incoherence leads the Fano resonance 
suppressed, and its asymmetric feature becomes consid- 
erably weak. First of all, shift of the peak position in the 
Fano resonance turns out not to occur in the multichan- 
nel Kondo effect. These features are argued to be quite 
general, distinguishing the non-Fermi liquid phase from 
the Fermi-liquid state in STM. 

The plan of the present paper is as follows. In Sec. II 
we present our model on an STM setup and derivation for 
the tunneling current. We introduce the nonequilibrium 
NCA in Sec. III. The conductance profile is analyzed at 
zero temperature in Sec. IV. Finally, the conclusion is 
presented in Sec. V. 



II. TUNNELING CONDUCTANCE 
A. Model 

The system under consideration is shown schematically 
in Fig. [1] It consists of a multichannel Kondo impu- 
rity host and a single- channel STM tip, placed directly 
above the host surface. The multichannel Kondo impu- 



rity host is modelled by a multichannel Anderson model 
in the slave-boson representation, which explicitly sepa- 
rates spin and channel excitations.— The STM tip cou- 
ples separately to the impurity and to the local conduc- 
tion electrons of the host. 

The Hamiltonian of the system takes the form 

H = -ffhost + -fftip + -fftunn (1) 

H host = ^ £ k4 CTT C k(JT + Sf Y ft fa 
k,<r,r tr 

+V c J2tihc kaT +h.c, (2) 

k,(j,r 

Hti P = ^(EK-eV^a^, (3) 

k,<7 

fftunn = ^ V a (Rt)fvbfa ka + h.C. 

k,(T,T 

+ *c(k,R||)4 CTT Sa+h.c., (4) 

k,p,(T,r 

where c£ {a ka ) creates host (tip) conduction electrons 
with wave vector k, spin <r, and channel r. The spin 
degeneracy is N and the number of channels is M. In 
the slave boson representation the impurity creation op- 
erator is given by f^bf, where f a is a fermion operator 
and bf is a boson operator. The fermion, /J, transforms 
according to SU(N) and creates a local spin excitation, 
whereas the boson bf transforms according to the conju- 
gate representation of SU(M) and annihilates the chan- 
nel quantum number of the "vacuum" state produced 
by destroying a conduction electron . 16 ' 17 Completeness 
of local states at the impurity site is represented by the 
constraint J^a /1/<7 + Et b\b f — 1, implemented as usual 
by introducing a Lagrange multiplier A. £k is the band 
dispersion of conduction electrons in the host, and Sf is 
the energy level of the Kondo impurity or the adatom 
at the host surface. V c is the hybridization parameter of 
the impurity and conduction electrons in the host. For 
simplicity, V c is assumed to be constant. E k is the band 
dispersion of tip-conduction electrons, and eV is an ap- 
plied voltage bias between the tip and host, which causes 
a weak electron current to flow between them. We set 
the chemical potential of host conduction electrons as 
our reference energy. 

Couplings between the impurity and tip-conduction 
electrons are represented by the hybridization parame- 
ter 14 (Rt), where Ri is the tip position. It decays with 
a tip-to-impurity separation 7 - 

F a (Rt) « Ke-^l, 

where k is an effective decay constant evaluated for states 
at the Fermi level of the tip. In this paper we consider 
\Rt\ <C 1/k, thus V^(Rt) is modelled as a constant. Cou- 
plings between the tip and host conduction electrons are 
represented by f c (k, Rii), where Rm is a parallel distance 
between the tip and impurity. For plane waves of con- 
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FIG. 1: Scanning tunneling microscope (STM) device with 
a tip placed closely to a Kondo impurity on the surface of a 
normal metal (host). In the host the impurity is hybridized 
with the metal conduction band through coupling V c . The 
STM tip couples to the impurity via hopping V a and to the 
local conduction electrons of the host via hopping t c . 



duction electrons, we have 

t c (k,R||) =i c e- ikR ii. 

When the tip is placed directly on top of the impurity, 
we take t(k, R||) = t c . 

The host Hamiltonian of Eq. ([2]) is the multichannel 
Anderson model. In the case of M = N the impurity 
is completely screened to form the Kondo singlet, re- 
sulting in the local Fermi liquid. When the number of 
channels is larger than the spin degeneracy (M > N), 
the impurity is overcompensated to give rise to a non- 
Fermi liquid fixed point, which exhibits universal power- 
law scaling.— _ — The universal scaling property lies at 
the heart of quantum critical phenomena in a number 
of materials. In this respect the present study can be 
said to probe non-Fermi liquid physics with STM. Un- 
fortunately, experimental realization of the multichannel 
Anderson model is still problematic. Recently, the two- 
channel Kondo effect was realized artificially in quantum 
dotSfii If the STM tip is applied to one of the leads, such 
non- Fermi liquid physics would be observed. 



B. Tunneling current 

The electron current flowing between the tip and 
host is calculated within the Keldysh nonequilibrium 
formalism i 18 ! 19 The current from the tip to the host is 
given by the time evolution of the occupation number 
for the electrons in the tip 

■Jt-fcW = e(^) - -j([H(t),N t (t)]), (5) 
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where N t = Eker a kcr a kcr- ^ne can express the current 
via nonequilibrium Green functions 

•w*) = J E KG&r^&t) 

k,cr,r 

+ \ E ^U^('- ( )+ k - ( fi ) 

k,p,cr,r 

where 

Gt T , a ^{t,t') = i{at a {t')d ar {t)), (7) 

G^a^frt) = »(oL(Ocp«rr(*)>, («) 

are lesser Green functions. Here we use the notation 

In the steady state nonequilibrium Green functions de- 
pend only on t— t', and the Fourier transformation results 
in 



and / a ( c )(w) is the Fermi-Dirac distribution function for 
tip (host) conduction electrons. T and the coefficients 
Qr(i) are defined as 



T 
Qr 
Qi 



4", 



(1 + M 7 ) 2 ' 
1 - M 7 



(I + M7) 3 
4 1 

(1 + m 7 ) 3 Y s 

4(1 -M7) 1 

(1 + m 7 ) 4 r~ 

r c ((l + M 7 )(7(M-l) + l)- 7 )) + 

(r a + 7 r c )(r c - 7 r a ) 

4(M- 1)7 1 



.(r a + 7 r c )(r c + M 7 r ) 
(r a - 7 r c )(M 7 r a + 



JE / ''-'Ov.,, Tr r u) 



(1 + m 7 ) 4 r. 



■(r c + 7r a )(2r c + (A/ 7 -i)r Q ), 



+ ^/<L)t;c» + ^. (9) where r s = r Q + r c ( 7 (Af-i) + i). r 

J is the coupling strength between tb 



where G daTaa (uj), Gf aT aa (uj) are the Fourier transfer- 

mations of Ek^^t.f), E k , P ^,ak.(M'), re- 
spectively. These nonequilibrium Green functions can be 
expressed via the impurity Green function based on the 
equation of motion method. Detailed calculations are 
presented in Appendix A. 

In a similar way we can calculate a current flowing 
from the host to the tip 



\H(t),N c (t)]) (10) 



t u\ / dNc 



with N c = Eko-r c kcrr c k(TT' where its detailed expression 
is given by Eq. (|A13I) in Appendix A. 

Calling Jh^t — —Jt^h in the steady state, the steady 
current can be rewritten in the form 



J = yJt^h - (1 - y)Jh^t, 



(11) 



where y is an arbitrary number. We choose y such 
that the term associated with the lesser Green function 
G drTT daT {uj) vanishes in the current formula from Eq. 
(|A12|) and Eq. (|A13I) . As a result, the electron current 
flowing between the tip and the host reads 

■/ = rV f ^ T tr (w) - /cH] , (12) 



<JT 



where 



T tr (w) - T + Q R ReGl TidaT (u) + QjImG^^), 

(13) 



a(c) = \y a (c)\ 2 ^Pa(c) 

is the coupling strength between the impurity and tip 
(host) conduction electrons, and 7 = ir 2 \t c \ 2 p a p c is a 
measure of the strength for the direct tunneling of con- 
duction electrons between the tip and the host. p a M is 
the density of states for noninteracting tip (host) conduc- 
tion electrons at the Fermi level. 

The current formula in Eqs. (|12p and (fl"3"|) can be 
viewed as a generalization of the Landauer-Buttiker for- 
mula to the STM case ) 24 ' 25 where T tr (o;) is the transmis- 
sion probability of the electron tunneling. The first term 
To of the transmission probability is the direct tunnel- 
ing between the tip and host, whereas the rest describe 
both indirect tunneling of conduction electrons through 
the impurity and interference between the two ways of 
electron tunneling. For the single-channel case (M = 1) 
the current formula in Eq. (fl2|) is reduced to the well- 
known formula.— In this case one may expect 7 <C 1 due 
to weakness of the tip coupling, hence Qr never van- 
ishes. For the multichannel case 7 = 1/M may happen 
when the channel number M is large. In this special case 
Qr = and interference contributions to the tunneling 
current vanish. When 7 = 0, i.e., there is no direct tun- 
neling between the tip and host, only the last term of the 
transmission probability in Eq. (|13l) appears, associated 
with the indirect tunneling of electrons between the tip 
and host through the impurity. This contribution is pro- 
portional to the density of states (DOS) of the impurity. 
In the Kondo regime the impurity is screened by conduc- 
tion electrons of the host, and this many-body effect must 
reflect in the DOS of the impurity, hence also in the tun- 
neling current. In general, the transmission probability 
T tr (a;) is a superposition of the continuous direct tun- 
neling, indirect tunneling, and their interferences, giving 
rise to the Fano resonance. 4 
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Discussion 



The linear conductance is given by 



G(w) = 



dJ(u) 



deV 



eV=0 



(14) 



at zero temperature. In general, the presence of the 
STM tip could affect physical properties of the host with 
the Kondo impurity We will discuss this point in Sub- 
sec. IIV CI However, if couplings of the tip to the host 
and impurity are weak, influences of the tip on the host 
and impurity are negligible, where the impurity Green 
function can be evaluated without couplings of the tip. 

The impurity Green function can be written in the 
form of the Dyson equation 



g: 



i 



ar.dar 



(15) 



where is the retarded self-energy of the impurity 

Green function. Considering real and imaginary parts 
of the retarded self-energy, E(w) = S_r(w) — iSj(w), we 
obtain the conductance profile in the form of 



T tr (cj) = T 



n 2 (uj) + 1 : 



(16) 



where 



q(w) = 



r c + E/(^) ' 

(1 + M 7 ) 2 Q R 



8"/ 



p(uj) = 1 - q 2 (uj) - 



(1 



- M 7 ) 2 



4 7 r c + S/( w )' 



i 1 i 1 i 1 i 1 i 

q=0, p=0.5 

■ q=0.5,p=0.5 " 

q=l, p =0.5 
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FIG. 2: Conductance profile in the noninteraction case for 
various parameters q and p as indicated in the figure. 



In the large frequency limit S_r(w) — > const, and 
—> result, thus we have T tI (u) — > T , nothing 
but the background profile of the conductance given by 
the direct tunneling probability between the tip and host. 

The conductance profile in Eq. ([To]) can be viewed as 
a generalized Fano form. In the noninteraction case of 
£(u;) = we have Q(oj) — (lu — e c ;)/r c , and p(lo), q(u>) 
are constants. In this case the Fano resonance results 
from the interference effect of a Lorentzian line-shape of 
a discrete level with a flat continuous background. The 
quantity q is the so called asymmetry parameter of the 
Fano line-shape, whereas the quantity p shifts positions 
of the maximum and minimum in the Fano line-shape. 
In Fig. [5] we present the conductance profile for various 
parameters of q and p. It shows the Fano-resonance line- 
shape like the Lorentzian one, the width of which is of 
order of r c . For q = the profile line-shape is symmetric, 
and its asymmetry becomes obvious as q increases. The 
parameter p not only affects the maximum and minimum 
positions of the Fano line-shape, but also makes the line- 
shape asymmetry clearer. One can expect that when 
interactions are included, the conductance profile will be 
significantly modified by the impurity self-energy as well 
as the Fano resonance. 



III. NONEQUILIBRIUM NONCROSSING 
APPROXIMATION 



The impurity Green function is evaluated within the 
nonequilibrium NCA, derived in a similar way as the 
equilibrium caseJ^ We start from an Anderson model 
in the slave-boson representation 

+ f dt^fi(t)(idt-e d -iX)U(t) 
+ J dt^b\{t){id t -i\)b?(t) + J dtiX 

dtV c f t(t)b f (t)c aT {t) + h.c, (17) 



where G Q caT caT (t, t 1 ) is the nonequilibrium single-site 
noninteracting Green function for host conduction elec- 
trons and time integration is performed along the 
Keldysh time contour. 

Integrating over conduction electron fields and intro- 
ducing two bi-local fields £/(£, t') and £&(£,£') conjugate 

to E<rfKt)U(t') and Er4(*)M*0» respectively, the 
quartic term in the effective action can be decoupled as 
follows 

S* cff = [dt^f^t){iat-e d -iX)U{t) 
+ J dtJ2 b Ht)(idt - iX)bf(t) + J 



dtiX 
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J J dtdt'Zfit' ,t)^2bt(t')br(t) 

J J dtdt'^it^Y^fi^f^t 1 ) 

dtdt'E b (t, t')D^{t, 0, (18) 



where D (t,t') — \ V c \ 2 G Q caT caT (t, t') is the hybridization 
function. 

The noncquilibrium NCA is the saddle-point approx- 
imation of the effective Keldysh action for the bi-local 
fields T,f(t,t') and Introducing nonequilib- 

rium fermionic and bosonic Green functions as F(t, t') = 
-i{T c f a (t)fl(t')) and B(t,t') = i(T c b f (t)4(f)), we find 
the saddle-point equations for the fermionic and bosonic 
self-energies 

£/(*,*') = iMD (t,t')B(t,t'), (19) 
E 6 (t,0 = -iND (t',t)F(t,t'), (20) 

where the bi-local fields play the role of self-energies of 
the fermionic and bosonic Green functions in the saddle- 
point approximation, given by 

F-\t,t') = 5{t,t'){id t -s d -iX)-E f (t,t'), (21) 
B-\t,t') = 5(t,t')(id t -i\)-Z b (t,t'). (22) 

Variation of the effective action with respect to the La- 
grange multiplier A gives rise to the constraint equation 



N(ftf 9 ) + M(bUr) = I- 



(23) 



Using the Langreth's rule of analytical continuation on 
the real time axis , 22 ' 23 the self-energy equations (fT9|) -([20 |) 
are 

Ef(M') = iM[(D R (t,t')+D<(t,t'))B R (t,t') 

+D R (t,t')B<(t,t')], (24) 

Ef(M') - -iN[D<(t',t)F R (t,t') 

+D A (t',t)F<(t,t')], (25) 

T,f(t,t') = iMD^(t,t')B < (t,t'), (26) 

£<(M') = -iND>(t',t)F<(t,t'). (27) 

In the steady state the Fourier transformation for the 
Green functions and their self-energies results in the fol- 
lowing nonequilibrium NCA equations 

£?M = MT C J ^-B<(e) 

+MT c f^f c (e-oj)B R (e), (28) 
SfH = NT C J ^-F<(e) 

+NT c J^-f c (s-oj)F R (s), (29) 



-N 



de 



f c {u-e)F<{e), 



(30) 
(31) 



where we have used explicit expressions for the hybridiza- 
tion function 

D R {u>) = |V c | 2 G°^ Cff » = -iT c , (32) 

D<{u) = \V c \ 2 G°< T ^ T (u}) = 2iT c f c (u}). (33) 

Note that the first terms in Eqs. (|28p and (|2"9")l are just 
constants. They can be absorbed into the Lagrange mul- 
tiplier, using the constraint equation (|23[) and 



(34) 
(35) 



We also used the fact that the Lagrange multiplier takes 
a large value at the end of calculations j^2i 

The Dyson equations (|2"Tj) - (l2"2"T) can be also rewritten 
for the retarded and lesser Green functions based on the 
Langreth's rule of analytical continuation 

F *M = \ ^R? \ ' (36) 



u) — Ed — i\ — Y, r (uj) 



B r {uj) 



1 



(37) 



to — i\ — T, r (uj) ' 

F<(u) = F r (uj)J:<(uj)F a (lj), (38) 

B<(u) = B r (lj)J:<(lj)B a (u J ). (39) 

Finally, the impurity Green function can be calculated 
via the fermionic and bosonic Green functions 

G R ^ daT (u } )=i J ^.[F<(u>+e)B A (e)+F R (u+e)B<(e)]. 

(40) 

Inserting this impurity Green function into the transmis- 
sion coefficient, we find the conductance profile measured 
in STM. 

The present derivation of noncquilibrium NCA equa- 
tions can be viewed as the path integral version for the 
projection method ) 20 ' 21 completely equivalent with each 
other. In practice, such NCA equations are first solved 
for retarded Green functions, and lesser Green functions 
are found with the use of the retarded Green functions. 
In the next section we will perform this work. 



IV. FANO RESONANCE IN THE 
MULTI-CHANNEL KONDO EFFECT 

A. Zero temperature solution of the noncrossing 
approximation equations 

In the linear response regime the host-electron distri- 
bution function f c (uj) is given by the standard Fermi- 
Dirac distribution function. Then, the NCA equa- 
tions (|2"5)) - ([2T))) for the retarded Green functions resemble 
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the equilibrium NCA equations, solved exactly at zero 
temperature^^ Eqs. (|58" |) -P9" j) can be written as the 
following differential equations at zero temperature 



dSf(u) 


MT C 


duj 


7T 


d££(w) 


Ar c 


du 


7T 



(41) 
(42) 



with the boundary condition Tif(-D) = E R (-D) = 0, 
where D is the band cutoff. 

Solving these NCA equations, one can find 



Ea-U) T, 



K 



t k ' r c 

where the scaling function h a [x] is given by 



h a [x,c] 



dy 



W[y° 



1 



(43) 



(44) 



W[y a e™°y] ' 

W[x] is the Lambert W function, defined as 28 

x = W[x] exp(W[a;]). 

T K = D[MT c /ttD] m / n exp[we d /NT c ] is identified with 
the Kondo energy scale, below which the multichannel 
Kondo effect arises. E the ground state energy of the 
impurity, below which spectral densities of the fermionic 
and bosonic fields vanish at zero temperature ! 26 ' 27 There- 
fore, F r (uj), B R (u>) are real functions below E . As 
shown in this expression, there are two energy scales T c 
and Tk for the NCA solution. Detailed derivation can 
be found in Appendix B. 

Equation (j43|) shows that the NCA solution obeys the 
universal scaling form. Although the scaling function in 
Eq. (|44|) should be computed numerically, its asymptotes 
in limits i«l and x S> 1 can be found analytically. For 
x <C 1, W[x] — x — x 2 , thus we obtain the asymptote 



h a [x,c] = [{a + l)^] 1 /^ 1 ) 



f - 



a + 2 



+ \(a + l)x} a/{a+1) 

2a + 1 



[(a+l)a:] 1/(a+1) 
(45) 



The leading term of the scaling function h a [x,c] shows 
the power scaling with an exponent I /(a + 1), imply- 
ing that both the retarded fermionic and bosonic Green 
functions exhibit power-law physics near the threshold 
energy Eq. This corresponds to the overcompensated 
regime, where the impurity spin is over-screened by mul- 
tichannel conduction electrons^ - — In the opposite limit 
x 3> 1, we obtain h a [x, c] = x, leading the fermionic 
Green function behave like l/w. This corresponds to the 
free moment regime, where the impurity spin is weakly 
bound to screening clouds. 

In Fig. [3] we plot the scaling function h a [x,c\. This 
NCA solution [Eqs. (|B6l)- (|43|) ] may be viewed as the 
complete solution of the NCA equations at zero tem- 
perature for all energy scales below the high-energy cut- 
off. It shows that the scaling function h a [x,c] crosses 




r 5 ( ^ 
logix) 



FIG. 3: (Color online) Scaling function h a [x, c] with a = 
0.5 and various parameters c. The two asymptotes [(a + 
l)x] 1 ^ a+1 * > and x are shown as the solid lines. 



from the power scaling regime to the linear behavior as 
x increases from zero. The scaling function h a [x,c*\ 
obeys the power law up to a; ~ 1, where c* is iden- 
tified with a "crossover" value. On the other hand, 
the power law of the scaling function is valid only for 
x <C 1 when c <§; c*. This property implies that there 
is a characteristic value of (Tk/F c )*, leading the power 
scaling to persist until energies comparable to Tk- In 
the conventional case of T^/r c <§; 1 the power scaling 
holds only for energies much below Tk- In the scaling 
regime F r {uj) - (|w - E \/T K ) M/{N+M) and B R {u) - 
f> - E Q \/T K ) N/(N+M) result, obtained previously^^ 
Although such power-law scaling fails to describe Fermi 
liquid in T < Tk for the single-channel case, it is the un- 
derlying physics of the overcompensated screening impu- 
rity in the multichannel Ccise, cis shown by the conformal 
field theoryJ^ 

The lesser self-energies obey the following differential 
equations 



dS<(w) 


MT C 


duj 




dZ<(u>) 




du> 


7T 



C F<H 



(46) 
(47) 



at zero temperature. In the scaling regime the lesser 
Green functions also display the same power scaling as 
the retarded ones, given b y 26 ' 27 



iA- 



1 



B < {uj) = -iA 



1 



Y b (cj) 



(48) 
(49) 



with A = 2tt/(N + M). See Appendix B. Based on 
Eq. (HD with Eq. (|B6]) and Eqs. @HJ), (gSJ}, we find the 
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final expression for the impurity Green function from 
Eq. P0|) in the scaling regime 



dar^dcrr 
1 



(0) + AG dtTT>daT (oj), 



(50) 



■ cot 



N + M 

r.dar 



N + MT C 
/ ttM \ 
\N + MJ 
in 

— i 



1 

M 



N + M 



N 



7T 



(N + M)T C 
M / 2N M 



N + M 
■ sin 



ttM 



2N + M 
ttT k 



N 
B 



M' AT- 
iV 



)( 



2M 



(N + 2M)T C \N + M N + M 



N + M. 
N + M u) 
M T~ K 
N + M uj 
M T) 



K 



where B(x,y) is the beta function, 29 and nj = 

The impurity Green function exhibits two power scal- 
ings with N/(M + N) and M/(N + M). In the overcom- 
pensation case (M > N) the power N/(N + M) scaling 
is dominant for uj < Tnca, where Tnca is the crossover 
energy when the dominant scaling behavior of the impu- 
rity Green function crosses from one power to another, 
given by 



Tnca = Tk 



M 



M (N + 2M)T C B {WTm> Tf+Ta 



(■ 



2N 



M 



N + M 
ttM 



(2N + M)ttT k B f 



N 



2AI 

N+M> N+M 



X COS 



N + M 



N + M 
M-N 



(51) 



B. 



Conductance profile within the noncrossing 
approximation 
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FIG. 4: Conductance profile calculated within the NCA (the 
solid lines), the marginal Fermi-liquid theory (the dashed 
lines), and the Fermi-liquid theory (the doted lines) for vari- 
ous values of 7. Other parameters are M = 6, N = 2, r c = 1, 
T a = 0.01, T K = 0.01, and n f = 0.8. 



The impurity Green function Eq. (|50[) is nonanalytic 
at uj = 0, exhibiting an asymmetric and sharp cusp with 
power-law scaling around uj — 0. These non-Fermi liq- 
uid features are reflected on the conductance profile in 
Eq. (IT6|) . For comparison, we also calculate the conduc- 
tance profile in the Fermi-liquid phase, given by the fol- 
lowing self-energy^ 

£ FI » = - r c [£-Hl(£) 2 ]. (52) 

In Fig. [4] we plot conductance profiles in the sym- 
metric case, i.e., = — ReE(0), within the NCA and 
Fermi- liquid theory for various values of 7. It shows that 
the conductance profile of the overcompensation multi- 
channel Kondo model shows a sharp cusp with power-law 
scaling at uj = Q, as expected. This feature is completely 
distinguished from the Fermi- liquid theory result, where 
the conductance profile exhibits the narrow Fano-Kondo 
resonance, the width of which is of order of Tk- Note that 
the frequency in the x-axis is scaled with Tr in Fig. 21 
which is much smaller than the energy scale T c in the 
noninteraction case (Fig. [2]). 



In the case of 7 — 1 /M the asymmetry parameter q(uj) 
vanishes, thus interference contributions to the conduc- 
tance disappear. As we can see in Fig. 01 the Fermi- liquid 
conductance profile is symmetric and exhibits only the 
Kondo resonance at zero frequency. On the other hand, 
the asymmetry parameter q(ui) is finite for 7 7^ 1/M, 
and the Fano resonance shifts the peak position away 
from zero, giving rise to an asymmetric feature in the 
conductance profile within the Fermi-liquid theory. Note 
that the asymmetry is more pronounced for larger 7. 

In the multichannel case the conductance profile still 
exhibits an asymmetric feature even when 7 = 1/M. 
However, this asymmetry is not due to the Fano reso- 
nance, but due to the non- Fermi liquid feature in the den- 
sity of states of the overscreened impurity. Although in- 
terference contributions to the conductance are finite for 
7 7^ they cannot shift the peak position away from 

zero in the multichannel conductance as in the Fermi 
liquid theory, the hallmark of non-Fermi liquid physics 
measured in STM. The conductance profile still exhibits 
the sharp cusp at uj = 0. This asymmetry is due to both 
interference contributions and non-Fermi liquid proper- 
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ties of the overscreened impurity. However, interference 
contributions are so weak that the Fano resonance is sup- 
pressed. This indicates dominance of the non- Fermi liq- 
uid overcompensation over the Fano resonance. In the 
multichannel case the conductance profile shows weak 
dependence on 7, thus the strength of the tip coupling 
is not essential to detect the non-Fermi liquid feature in 
the tunneling conductance, as far as it does not vanish. 



Physical origin for suppression of the Fano 
resonance 



In the Fermi liquid phase one can see Sl(w) ~ uj 
and q(co), p(u>) ~ const as w -4 in the tunneling 
conductance [Eq. (fTBT) ] from Eq. (|52|) . giving rise to 
the Fano resonance away from u> = 0. This originates 
from quantum coherence of impurity dynamics, which 
maintains Efl(w) ~ uj and Sj(cj) ~ uj 2 as uj — > 0. 
In the multichannel overcompensation scaling regime 
E R (u) ~ - \ U \ N /( N + M ) results when w _> 0, 

which leads tt(uj) ~ 1 and the asymmetry factor q(u>) ~ 
\u\- N /( N + M ), p(u) ~ \uj\~ n ^ n+m \ As a consequence 
we find T tl (w) ~ \ W \- 2N /( N + M ) in w ^ . The conduc- 
tance profile always exhibits its peak at to = that indi- 
cates the peak position does not shift unlike the Fermi- 
liquid theory. This can be interpreted as the fact that 
quantum coherence of impurity dynamics is lost in the 
multichannel scaling regime, hence suppressing the Fano 
resonance. 

For comparison we also calculate the conductance pro- 
file in a marginal Fermi-liquid phase, given by the follow- 
ing self-energy ansatz^L 



D. Effect of the tip- host and the tip-impurity 
coupling on the M-channel Kondo impurity system 



JMFL 



(w) = r c 



W , U 7T \oj\ 

log I 

T K & T k 2T k 



(53) 



The marginal Fermi-liquid phase can be considered as 
an intermediate case which lies between the Fermi-liquid 
phase and multichannel overcompensation phase, where 
the impurity dynamics maintains weak quantum coher- 
ence. In Fig. [4] we also plot the conductance profile in 
the marginal Fermi-liquid phase for various values of 7. 
It shows that the conductance profile of the marginal 
Fermi-liquid theory exhibits a narrow resonance nearby 
uj = 0. It resembles to the conductance profile of the 
Fermi-liquid theory, except for the position of the Fano 
resonance which very closes to ui = 0. The weak quan- 
tum coherence in the marginal Fermi-liquid phase can 
maintain the Fano resonance, however, weakness of its 
quantum coherence cannot significantly shift the peak 
position of the Fano resonance away from zero in com- 
parison with the Fermi-liquid phase. This feature implies 
the important role of quantum coherence of impurity dy- 
namics in the mechanism of Fano resonance. When the 
impurity dynamics becomes incoherent, the Fano reso- 
nance is suppressed. 



As mentioned in Subsec. fll C[ the impurity Green func- 
tion is evaluated without the tip coupling in the conduc- 
tance profile of Eq. (TIB)) . In general, the presence of the 
tip could affect physical properties of the host with the 
Kondo impurity. In this subsection we discuss the effects 
of the tip-host and the tip-impurity coupling on the M- 
channel Kondo impurity system of the host and impurity. 
Although such tip couplings break the SU(M) symmetry 
in principle, we argue that the overcompensation Kondo 
effect still occurs for weak tip-impurity couplings with 
finite tip-host couplings. 



1. Effect of the tip-host coupling 

First, we consider the effect of the tip-host coupling 
only, where the tip-impurity coupling is neglected. We 
take the following unitary transformation for conduction 
electron fields Cko-r and bosonic fields b T 



Ck 



— ^ T_J rr ' Ckcrr' j 
t' 



(54) 
(55) 



where the M x M unitary matrix U is chosen to satisfy 
Ckcri = St c kcrr/V^I for diagonalization of the tip to 
host coupling term. We rewrite the starting Hamiltonian 
Eq. (fT]) in the above transformed basis and integrate over 
tip conduction-electron fields. Then, we find 

Scs = J dT^2cl aT (r)(d T -e k )c kt7T (T) 

k(TT 

+ E f^ d r + X ~ F ~f)U T ) + E ^)(5t + X)b f (r) 

<T T 

+ KE^( T )M T )WT)+h.C. 
kcrr 

+ f dTdr'M\t c \ 2 J2ct al ( T )9a(r~T')d k , t7l (T'), 

J kk'cr 

where g a (r — t') is the local Green function for tip elec- 
trons. 

It is clear that the last term in the above effective 
action breaks the SU{M) symmetry. Effectively, the 
conduction channel Ck<ri has an additional contribution 
~ M\t c \ 2 p a to its normal dispersion, while other con- 
duction channels have not. Integrating over conduction 
electrons, we obtain 

Scff = f dr^2ft(r)(d T + A - e f )Mr) (56) 
+ J2U(r)(d r + A)6t(t) + JdTdr'\V c \ 2 
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E^)^){E^(t-tO}&I(tO/.(t') 

a kk' 

+ J drdr'\V c \ 2 £ /J(r)M^){E&'( T ^)} 

o-,-r#l kk' 

^OWA (57) 

where [^(r-r')]- 1 = -(5 T -£ k )5 fcfe - -M|i c | 2 5Q (r-r') 
is the inverse of the electron propagator for the channel 
t = 1 and [g% k t{r - r')] _1 = -(d T - Ek)5 kk , is that for 
other M — 1 channels. This SU(M) symmetry breaking 
gives rise to anisotropic hybridization couplings. In Ap- 
pendix C we prove that T { <t ] = |F c | 2 £ kk , l Im 3fefeV (0)1 < 
r c always happens. This indicates that the channel r = 1 
couples to the impurity weaker than (M — 1) rest chan- 
nels. As a consequence, the [M — 1) channel Kondo ef- 
fect would occur. For large M, there is no difference 
in physics of the overcompensated Kondo effect between 
(M — 1) channels and M channels. 



2. Effect of the tip-impurity coupling 

Now we consider the effect of the tip-impurity cou- 
pling only, where the tip-host coupling is neglected. As 
performed in the previous subsubsection, we take again 
the unitary transformation in Eqs. (|54|) -(|55 j) for con- 
duction electron fields Cko-r and bosonic fields b T , but 
now the M x M unitary matrix U is chosen to satisfy 
bi = J2 T b T /VM for diagonalization of the tip to impu- 
rity coupling term. Then, we find the following effective 
action 

S cS = [ drJ2ft(r)(dr + X-e f )U(r) 

J a 

+ J2bHr)(d T + X)b f (r) + f drdr'Y,fl{r)Hr) 

t •'a 

(\V c \ 2 9c (r t') + M\V a \ 2 g a (r t'))1\{t') U{r') 

+ J drdr'\V c \ 2 £ fl{T)b T {T)g c (T~T>)b\{T>)UT'\ 

where g c (T) = ^ k k' 9kk' ( T )- This effective action shows 
that the tip-impurity coupling breaks the SU (M) sym- 
metry. The effective hybridization coupling of the ef- 
fective channel r = I is T c + MT a , and it is always 
larger than the effective hybridization coupling T c of the 
{M — 1) rest channels. As a consequence, the single- 
channel Kondo effect would occur. In this case the 
system restores the Fermi-liquid behaviors. The NCA 
at zero temperature produces spurious non-Fermi liquid 
features ) 26 ' 27 and it is not adequate to describe the physi- 
cal properties of the system. In the calculation of the con- 
ductance profile the Fermi-liquid self energy in Eq. (|52|) 
must be used, and the Fano-Kondo resonance could ap- 
pear. 



3. Effect of both the tip-host and tip-impurity couplings 

The two previous subsubsections show that the tip- 
host coupling allows the overcompensation Kondo effect, 
while the tip-impurity coupling drives the system away 
from the critical regime. Therefore, we are mainly inter- 
ested in the case of a finite tip-host coupling and small 
tip-impurity coupling. As the case of finite tip-host cou- 
plings, we take the unitary transformation for conduction 
electron fields c ko -r and bosonic fields b T in Eqs. (j54|) - (j55|) 
with Cko-i = X^r Cko-r/v^f- Proceeding in a similar way 
as the previous subsubsections, we find the effective ac- 
tion 

S cS = f dTj2rt(r)(dr + \-e f )UT) 

+ J2bHr)(d T + \)b f (T) + J2 [drdr'\V c \ 2 

ft(r)b T (r)g c (r - t')1\{t')U{t') + W drdr' 

{\V c \ 2 fUr)b 1 (r)g^Hr-rrb\(r')Ur , ) + 
P a (T)b (r)[M\V a \ 2 g a (r - r') + M 2 \V a \ 2 \t c \ 2 J dr x dr 2 
9a(r - rO^Cn - T 2 )g a (r 2 - r')] &J(r')/ CT (r') + 
MV c V:t c fUr)bi(r) J d n g c ^(r - n)g a (n - r') 

bl(T')f a (r') + MV:V a t:P a (r)b (T) J dng a (r - n) 

g cW (n-r')b[(r')Ur')}, (58) 

where g°W( r ) = E k ,k'4'V)> and b ° = E T ^/VM. 
In general, &o is not necessarily equal to b\. We can 
rewrite bo = aibi + a 2 b, where b is a linear combination 
of b 2 , ... , 6m- The coefficients oc\ and a 2 satisfy |ai| 2 + 
|a 2 | 2 -l. 

In the wide band limit of a constant density of states 
for conduction electrons in both the tip and host, the 
effective hybridization coupling terms of M channels in 
Eq. ([55]) can be rewritten schematically as follows 

/^(rw + K| 2 ri ) + a \T ac + ai r* ac )blf a 
+/]6| a2 | 2 ri°)St/ CT + Jtu 

+rtba 2 T*J\f a + E ffr-T&rU, (59) 

where 

r(°) = M|y a | 2 |Im 5a (0)|+M 2 |F a | 2 |t c | 2 |Im[ 5a (0)V (1) (0)]|, 
T ac = MV c V:t c \Im[ 9a (0)g c ^(0)]\. 

One can show that T { c 1] = r c /(l+M 7 ), ri 0) = Mr a /(1+ 
Mi), and |r ac | 2 = M 2 r c r a 7/(1 + M7) 2 . The hybridiza- 
tion coupling of (M — I) channels T c is expected to dom- 
inate over those of the channels b± and b if the following 
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conditions are satisfied 



Acknowledgments 



r c > r«+if) + 2|r ac |, 
r(D > rl°). 



(60) 
(61) 



These conditions are equivalent to T c > MT a and jT c > 
(3 + \/8)r a . When the hybridization couplings of (M — 1) 
channels are relevant, the (M — 1) channel Kondo effect 
would occur. For large M, there is again no difference 
in physics of the overcompensated Kondo effect between 
(M — 1) channels and M channels. However, we admit 
that this expectation does not have a firm ground because 
there exists complex mixing between channels. In our 
opinion this problem should deserve further investigation. 

Finally, we would like to mention that our demonstra- 
tion is expected to apply to the quantum dot realization 
for the multi-channel Kondo effect. If the STM tip lies in 
the lead, we would be in weak T a . Then, the suppression 
of the Fano resonance is expected to observe. 



V. CONCLUSION 

In the present paper we predicted an interesting feature 
of non-Fermi liquid physics for the multichannel Kondo 
model, based on the STM experimental setting (Fig. [TJ) . 
Non-Fermi liquid physics often occurs at the quantum 
phase transition, and presents challenges in both theo- 
retical and experimental aspects, for instance, quantum 
criticality in heavy fermion materials*^ Heavy fermion 
materials can be modelled by the Anderson lattice model. 
One heavy fermion quantum critical point in the An- 
derson lattice model was argued to be captured in the 
so called dynamical mean-field theoretical framework, 
more concretely, the two impurity Anderson model with 
self-consistency, expected to result in essentially similar 
physics with the multichannel impurity model.™ As the 
first step to understand non-Fermi liquid physics within 
the STM detection, we employed the multichannel An- 
derson model for one source of the non- Fermi liquid state. 

We derived the Landauer-Biittiker formula for the tun- 
neling current from the STM tip to the multichannel im- 
purity host based on the Keldysh nonequilibrium formal- 
ism, where the tunneling current is given by only the 
impurity Green function. Employing the nonequilibrium 
NCA, we showed that the impurity Green function ex- 
hibits universal power-law scaling at low energies. As 
a consequence, the tunneling conductance turns out to 
exhibit weak asymmetry but rather sharp cusp at zero 
energy, resulting from the power-law scaling of the impu- 
rity Green function. The conventional Fano resonance in 
Fermi liquids was shown to be suppressed. The main pre- 
diction of our study is that the peak position in the Fano- 
Kondo resonance does not shift, even increasing the tip 
coupling constant, clearly distinguished from the Fermi 
liquid theory. Quantum coherence of the impurity dy- 
namics turns out to play an important role in the Fano 
mechanism. 
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Appendix A: Derivation of tunneling current 

In Appendix A we derive the tunneling current for- 
mula in Eqs. (fT2"|) and (IT3"1) . The lesser Green functions in 
Eq. (|9]) are given by off diagonal components of nonequi- 
librium Green functions, defined on the Keldysh time 
contour which runs on the time axis from — oo to oo, and 
goes back to — oo ) 18 ' 19 

G c dtTT , a Ut,t') = -i(T c d aT (t)al a (t')), (Al) 
^ pCTT , akCT (M') = -*<T cCpCTT (i)aL(0), (A2) 

where T c is the time ordering operator on the Keldysh 
time contour. Differentiating the nonequilibrium Green 
functions with respect to t' and resorting to the Heisen- 
berg equation of motion 



dt' h ] 



[H(t>),al(t>)]. 



we obtain the following equations 
GdaT^ka^it') = J ' dti V a G c daTdlJT {t,ti)g c aV .{ti,t') 

+ J2j dt ^ tcG daT ^ T ,(tM)9 C a*{tiX), (A3) 

GWp^foO - 5 TT , J dt x V c G d!TrMT {tM)9 c cp {tut') 

+ £^,aka(*.*l)<&,(*l,*'), ( A 4) 

k J 

where (^(t, <?ck(^') are ^ ne nonequilibrium Green 
functions for the isolated noninteracting conduction elec- 
trons in the tip and host, respectively. 

Inserting Eq. (|A4[) into Eq. (| A3[) ■ we obtain 



G daTiaa {t,t') — Jdti V a G c daTfdaT {t,tx)gl{tx,t') 

+ [dhjdh t c V c G daT dcrT (t, tl)</c(*l(*2)£f£(*2> *0 

IdtJdh M\t c \ 2 G daTMa {t, h)£(ti, t 2 )g c a (t 2 , tO,(A5) 



where g c a{c) (t,t') = £ k ffo(c)k(*> Usin § the Langreth's 
rule of analytical continuation on the real time axis ; 22 i 23 
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we can obtain the retarded (advanced) and lesser (or 
greater) Green functions from Eq. (|A5[) 



G 



R/A 



V G R/A * n R / A 



yR/A 

7 dar.dar 



*g? /A *g« /A 



(A6) 



a (G R aT 



dar * 9a 



Gdar,daT * 9a 



+t c V c (G R T *g R *g<+G, 



dar ..< 



■ 9c * 9a 



-M 



^ G d*T,d<?T *9c*9a 

\tc\ 2 (G% XT>ac *g**g„ ( T ;,; T7 . (((T i u, i ,/„• 



*5a 



(A7) 



where the superscripts, J?, A, <, denote the retarded, ad- 
vanced, lesser Green functions, respectively. For simpli- 
fying to write equations, we use the so called * notation, 
defined as A(t) * B{t') = f dtiA(t,h)B(ti,i/). 

In the steady state Eqs. (|A6|) - (|A7|) are easily solved, 
making the Fourier transformation. As a result, we find 



G 



dar.acr 



(t c V c + M\t c \ 2 )(g R {u)gf{u>)+gf{u)g A {u)) 



g: 



\d(TT ' 



(co)Z A (u) V a g A {w) + t c V c g A (u J )g A (u J ) ,(A8) 



where Z R ' A {u) = 1/[1 - M\t c \ 2 g R/A {u)g R/A {u)\. 

Performing in a similar way, we can express the 
Green function G^ a(7 (t,t') = Ek, P ,r ^w, opCT 
with G daaa .{t,t') = £ T G d ^ aa (t,t'). One can verify 

G< i0 » = Mt c Z R (u>)[g R (u)g<(u) + g<(u>)g A (u;)] 
[l + M\t c \ 2 Z A (u;)g A (u>)g A (u>)} 
+G< CT , a »Z» [V; 9 R [u) + Mt c V:g R {u)g R {u)] 
+G A r7 , a Jo J )Z R (o J )[v; 9 <(u) + {g R {u)g<{u>) + 

g<(u;)g A (uj))(Mt c V: + M\t c \ 2 Z A (u) 

(V c *g A + Mt c V:g A (u)g A (u))} . (A9) 



Using Eqs. (|A8|) - (|A9|) . we can express the steady cur- 
rent in Eq. ([9J with the nonequilibrium Green functions 
of the impurity. For simplicity, we will consider a flat 
density of states for the tip and host conduction elec- 
trons in the wide band limit, given by 



R/A, \ 



Ti^Pa(c), (A10) 
< (c) (w) = 2nip a(c) f a{c) (u}), (All) 



where p a i c ) is the density of states for noninteracting 
tip (host) conduction electrons at the Fermi level, and 
/o( c )(w) is its Fermi-Dirac distribution function. From 



Eqs. ©, 



, and (IA9I) we obtain 



d(TT,d(7T 



2tt I 2 

(Lu)T acl +G R d(r »2T Qcl / c (u;) 



where 



-G R aT>d(TT (uj)T ac2 (f a (uj) ~ / e (w)) +h.c.}, (A12) 



4" 



To 

T ac l = ITTpa 



(1 + M 7 ) 2 ' 

l V a + iV c t c irp c \ 2 



Tacl = 2inp a 



(1 + M 7 ) 2 ' 
(1 - MjYYa ~ iVctcVPc? 
(1 + M 7 ) 3 



Here 7 = tt 2 | t c 1 2 p a /o c is a measure of the strength for the 
direct tunneling of conduction electrons between the tip 
and host. 

In a similar way we can find a current flowing from the 
host to the tip based on Eq. (fTDj) . 



+G. 



d<JT,d<JT 



(u)T cal + G R ^ d<rr (cj)2T cal f a (u) 



+ 



+G£ rjdCTT HT CQ2 (/ c H - /„( W )) + h.c.} 

2G r t ^ t (uj) [AT c f c (uj) + ATafaiu) 
AT(/ c ( W )-/ ( W ))]+h.cj, (A13) 



where 



1 



1 + M7' 



AT 



ir, 1 



1 



1 + M7 



2ir r 



1 



(1 + M7) 2 



M7 - 1 

' (1 + Af 7 ) 2 



T^ai, T" ca 2 are just T ac i, T ac2 , changing indices a -o- c 
while the hopping £ c is unchanged. 

One can notice that the first term of the current Jh^t 
in Eq. (|A13[) is just the current J t ^ h in Eq. (|A12j) if the 
tip and host are interchanged with each other. In the sin- 
gle channel case (M = 1) the second term of the current 
Jh^t in Eq. (|A13[) vanishes, thus the current formula sat- 
isfies the symmetry between the tip and host. However, 
in the multichannel case (M > 1) the symmetry is bro- 
ken, because host conduction electrons are multichannel, 
whereas tip-conduction electrons are single-channel. 
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Appendix B: Solution for NCA equations 

In Appendix B we derive the NCA solution Eq. (|43|) 
with Eq. dHJ) and Eqs. g5]), (gU) from Eqs. flU]), ^ and 
Eqs. P5|) . (|4"T)) . respectively. Introducing inverse Green 
functions 



one can rewrite Eqs. (|41 [) and (1421) as 

dy/(w) , Mr 



dY b (uj) 



7T 

ivrv 



(Bl) 
(B2) 



(B3) 
(B4) 



Then, we find the exact relation between Yf(u) and Y b (ui) 

(B5) 

where T K = D[MY c /nD] M / N exp[-Ke d /NY c ] is the 
Kondo energy scale. 

Solving Eq. (|B5[) . we obtain 



Vr*/ v \mt c t k 



(B6) 



where W[x] is the Lambert W function, given by x = 
exp(VF[a;]). 28 Inserting Y b in Eq. (jB6| into Eq. (|B3t . 
we find the solution given by Eq. (|4"3"f with Eq. (|4"4")l . 

Although the NCA equations (|3S ]l -(|iT ]l for lesser self- 
energies are identical to the differential equations (|4ip- 
(142 p of retarded self-energies, the Dyson equations (|38jl- 
(f3T))) for the lesser Green functions have a different struc- 
ture. Together with the Dyson equations (f38|) - (f39|) . 
Eqs. (|4^1) - (|4T|) can be rewritten as 



d[F<(o;)y/M] 
du 

d[B<{ui)Y b \uj)} 
dw 



MY, 



IT 

NY 



B<{u), 



for frequency below Eq. Then, we find Eqs. 
(l49l) as their solution. 



(B7) 
(B8) 
and 



Appendix C: Effective channel couplings 

In this Appendix we calculate the hybridization cou- 
pling of the channel r = 1 r c W for the effective action 
of Eq. (|S7]). 

We start from the following matrix identity 



(R + zU)- 1 = R 1 



1 + xS 



(CI) 



where = Ridij is a diagonal matrix, = 1 for all 
i, j, a; is a number, S = Yli^ii and = R^ RJ . 
Taking R = 8kk'{u ~ £fc) and x = — M \t c \ 2 g a (uj), we find 
the inversion of [ffj^Vtw)] -1 , 

Skk' , Af|i c | 2 5a (o;) 



i l 



w - £/c 1 - M |i c | 2 g a (w)3 c (w) w - e k u - e k > ' 

where g c (w) — X)k £ fe)- Performing the summation 

of k, k' we obtain 



kk' 



M|i c | 2 <? a (cj) 



M|t c | 2 5a (w) 5c (w) 



l-M\t c \ 2 g a (u>)g c (u))' 



(ff C ( w )) 2 
(C3) 



In the wide band limit of a constant density of states for 
conduction electrons in both the host and tip, Eq. (|C3[) 
results in 



r c(l) 



1 + M7 



(C4) 



It also shows T C W < Y c . 

We prove the identity Eq. (ICip . Consider 



(R + 1U)- 1 = ^(-l)"a;"(R _1 U)"R _1 

We introduce the following identity 

00 

^(-l) n a;"(R- 1 U)"R- 1 = 

n=l 

00 

-a; I" ^ (—l) n x n (R _1 U) n+1 R _1 



(C5) 



n=0 



(C6) 



Since the matrix R is diagonal, we can write its inverse as 
(R-%- = RtSij, where R t = I/R4. Calling (R-^jy = 
Ri, where elements in each raw are identical, we find 



[(R _1 U) n ] . = RiS n ~ 1 



(C7) 



with S = Y,i R i- Inserting Eq. JC7| into Eq. ([Co) , we 
obtain 

00 00 
^(-l)"a;™(R- 1 U) n R- 1 = -x\^2(-l) n x n S n F 

(C8) 



71 = 

X 



1 + xS 



where Fy = RiRj. Resorting to Eqs. (|C5|) and (|C8I) . 
finally, we reach Eq. (|C1|) . 
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